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Abstract 

In this paper we investigate the motion of small compact objects in non-vacuum spacetimes 
using methods from effective field theory in curved spacetime. Although a vacuum formulation 
is sufficient in many astrophysical contexts, there are applications such as the role of the self¬ 
force in enforcing cosmic-censorship in the context of the overcharging problem, which necessitate 
an extension into the non-vacuum regime. The defining feature of the self-force problem in non¬ 
vacuum spacetimes is the coupling between gravitational and non-gravitational field perturbations. 
The formulation of the self-force problem for non-vacuum spacetimes was recently provided in si¬ 
multaneous papers by Zimmerman and Poisson [1] and Linz, Friedmann, Wiseman [2]. Here we 
distinguish ourselves by working with the effective action rather than the field equations. The 
formalism utilizes the multi-index notation developed by Zimmerman and Poisson [1] to accom¬ 
modate the coupling between the different fields. Using dimensional regularization, we arrive at 
a finite expression for the local self-force expressed in terms of multi-index quantities evaluated 
in the background spacetime. We then apply the formalism to compute the coupled gravitational 
self-force in two explicit cases. First, we calculate the self-force on a massive particle possessing 
scalar charge and moving in a scalarvac spacetime. We then derive an expression for the self-force 
on an electrically charged, massive particle moving in an electrovac spacetime. In both cases, the 
force is expressed as a sum of local terms involving tensors defined in the background spacetime 
and evaluated at the current position of the particle, as well as tail integrals that depend on the 
past history of the particle. 


I. INTRODUCTION 

Solar-mass compact bodies spiralling into super-massive companion objects are promis¬ 
ing sources of low frequency gravitational waves for future gravitational wave detectors. 
The possibility of measuring these signals has fueled research efforts to construct models 
describing the motion of the small body which go beyond the test-particle approximation. 
The particle’s motion is influenced by the perturbation it sources, causing its orbit to de¬ 
viate from geodesic motion in the background vacuum geometry [3, 4]. The dissipative 
component of this “self-force” acting on the particle causes it to lose energy and angular 
momentum to gravitational waves leading to a decay in the orbital radius and an eventual 
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plunge into the larger companion body after a time Tinspirai ~ The gravitational 

self-force in vacuum spacetimes has developed into maturity over the last decade with no¬ 
table achievements being: the rigorous formulation in which the compact object is treated 
as an extended body instead of a point particle [5, 6], the numerical computation of the 
self-force in Schwarzschild spacetime [7] and its effect on the innermost stable circular orbit 
of a Schwarzschild black hole [8], the orbital evolution of a massive particle in Schwarzschild 
spacetime [9], the self-consistent evolution of a scalar particle around a non-rotating black 
hole [17], the computation of the scalar self-force in Kerr spacetime [10, 11], and the exten¬ 
sion to second-order in the mass ratio [12-16]. In addition, the self-force was also successfully 
compared with numerical relativity [19] and provided useful information for post-Newtonian 
£t parameters [20]. 

For many applications, the black hole sits in isolation of external fields and matter, in 
which case the vacuum formulation of the gravitational self-force is adequate. Other appli¬ 
cations, however, may require an extension to non-vacuum spacetimes. For instance, one 
might consider the motion of a satellite around a material body having some companion 
matter distributed over the region of the spacetime surrounding it. A non-vacuum formu¬ 
lation will also be required to elucidate the role played by the self-force in scenarios that 
aim to produce a counter-example to cosmic censorship by overcharging a near-extremal 
Reissner-Nordstrom black hole [21]. Hubeny [22] has shown that a nearly extremal charged 
black hole can absorb a particle with parameters living in an open three-parameter family of 
charge, mass, and energy such that the final configuration’s charge-to-mass ratio exceeds the 
extremal bound leading to a naked singularity. Jacobson and Sotiriou [23] showed that this 
phenomenon is not unique to the Reissner-Nordstrom spacetime when they revealed that a 
near-extremal Kerr black hole can also absorb a particle and be driven beyond the extremal 
bound. These works treated the particle as a test particle in the black-hole spacetime, and 
it was soon realized that self-force and radiative effects can play an important role in these 
overcharging and overspinning scenarios. In fact, Hubeny incorporated approximate self¬ 
force effects in her original analysis, Barausse, Cardoso, and Khanna [24] took into account 
the gravitational radiation emitted by the particle on its way to overspin a Kerr black hole, 
and work by Colleoni and Barack suggests how the full self-force may be incorporated into 
an overspinning condition [25]. 

Up until recently, the formulation of the gravitational self-force has been restricted to 
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bodies moving in vacuum spacetimes. In non-vacuum spacetimes, the gravitational pertur¬ 
bation induced by the small body generates a small deviation in the stress-energy tensor of 
the external held which generates coupled held equations. The coupling between the held 
equations makes the problem more difficult but not altogether intractable and progress has 
been made [1, 2, 27, 28]. The hrst full treatment of the coupled self-force on point particles 
in non-vacuum spacetimes was simultaneously given in Refs. [1, 2]. In [1], regular self-force 
equations of motion were derived via the Detweiler-Whiting axioms for Lorenz gauge per¬ 
turbations of both electrovac and scalarvac spacetimes. The treatment in [1] made use of a 
condensed index notation which led to held equations and expressions for the regular and 
singular solutions for arbitrary external held and point-like source content. The present pa¬ 
per will attempt whenever possible to carry on the same notational conventions as [1]. The 
simultaneous publication by Linz and collaborators [2] presented the self-force on a charged 
particle in an electrovac environment. They regularized the helds via a combined approach 
that used an iterative scheme based on the structure of the perturbed held equations, angu¬ 
lar averaging, and mass/charge renormalization. They also derived an explicit expression for 
the gradient of the singular held and showed that the coupling between the helds does not 
contribute to the “singular self-force”; the regularization parameters for the held gradients 
in the coupled problem are simply the sum of the decoupled parameters. 

The aim of the present work is to provide an additional formulation of the hrst-order self¬ 
force problem for non-vacuum spacetimes using the methods of ehective held theory (EFT). 
The EFT program is well suited to the binary problem because of its inherent separation 
of scales. EFT was hrst successfully applied to the problem of motion for slowly-moving, 
weakly-gravitating bodies in the post-Newtonian regime [31], but it has also proven itself as 
a contender in the self-force arena through the work of Hu, Galley, and collaborators [34-37]. 
The ehective held theory approach to world line dynamics has several attributes including a 
systematic power counting scheme via Feynman diagrams, dimensional regularization, and 
internal self-consistency (no ambiguity in dehning the world line). In the EMRI scenario, 
the EFT approach treats the small compact body as an ehective point particle, which takes 
the form of a skeleton world line with a set of multipole couplings. The ehective action 
preserves the symmetries of the full theory, which in the present case are general covariance, 
reparameterization invariance, and rotational invariance. This description is sensible at 
length scales asymptotically greater than the characteristic size of the object, but it breaks 
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down at the body scale just as a far-zone expansion of a radiation field breaks down when 
approaching the near-zone. Information about the internal structure of the small object is 
incorporated by a matching procedure that involves adjusting the coefficients parameterizing 
hnite size effects in the point particle theory to agree with the values obtained using a 
complete theory describing its internal structure. But, as it has been shown that hnite size 
effects such as tidal deformations are irrelevant for self-force computations up to fourth order 
in the mass ratio [35], such couplings will be ignored. 

The main results of this work are the equations of motion for a small body due to the 
coupled self-force in several non-vacuum scenarios. In all cases we hnd that the equation of 
motion can be decomposed into the schematic form 

ma- = FS + + F" „ 

where Fq is the force resulting from the gradient of the background potential, is the 
local contribution to the self-force which is built from background quantities evaluated on 
the world line, and is the non-local contribution to the self-force which takes the form of 
a time integral over the particle’s past history. Here we are not interested in examining the 
solutions of the field equations. Instead, they are “integrated out” from the problem in the 
process of obtaining an effective action governing only the world line. As a consequence of 
this, the time-dependent mass of the scalar particle in the scalarvac scenario automatically 
includes the effects of the regular self-held. The form of the equation of motion is unaffected 
by these mass corrections because they are suppressed by an additional power of the small 
mass ratio. The explicit form of the time-dependent mass due to self-force corrections can 
be found in [1]. 

We begin in Sec. II with a toy problem involving a multi-component scalar held in curved 
spacetime which is coupled to a world line with N scalar charges. This model captures the 
general features of the more complicated scalarvac scenario and ahords us the opportunity to 
introduce the EFT procedure for arriving at the ehective action via the generating functional. 
We also introduce several technical tools used in computing the hnite self-force such as the 
quasi-local coordinate expansions of various background items in the bulk and around the 
world line. 

In Sec. Ill, we introduce a general formulation which uses the multi-index notation of [1] 
to derive the ehective world line action for a set of coupled perturbations while leaving the 
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field arbitrary. In Sec. IV, we derive a causal equation of motion from the effective action 
in terms of retarded Green functions. Finally, in Secs. VI and VII, we apply the formalism 
to the scalarvac and electrovac spacetimes and derive the coupled self-force equations of 
motion. Since its publication, a small error in the local scalarvac self-force was found in the 
results of [1]. The equations of motion are presented here in their corrected form. However, 
because we do not derive an expression for the singular/regular fields explicitly here, we are 
unable to do a direct comparison with the results of [2]. 


II. MULTI-COMPONENT SCALAR FIELD 


A. Effective action and formal equation of motion 


Before launching into the complicated coupled problem involving the gravitational field, 
we choose to examine a simple toy problem which ignores the gravitational field altogether, 
but which captures the features of the more general case. The model we consider is a 
multi-component scalar held $^( 2 ^) of dimension N which couples to a massive world line 
possessing N scalar charges qa- The bulk action for the theory is given by 


^buik = -^ / + 2V($)) dV, (1) 

where is an iV x iV dimensional matrix and V(<h) is an arbitrary potential even in 

its argument. The world line action of the theory is a massive point charge interacting with 
the scalar held according to 




pp 



( 2 ) 


The motion of the scalar charge in the background spacetime is dictated by 


+ uV) (3) 

where we have introduced an ehective mass m := mo — depends on the 

scalar held, as well as an Einstein summation convention for “internal” indices, q^^A '■= 
X]a=i Qa^a- The presence of the point particle sources a scalar perturbation /a- The particle 
will generally also source gravitational perturbations, but we postpone their consideration 
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until later in the paper. The quadratic bulk action is given by 

5'^bulk = -^ I + 

+ + V^^^^fAfB)dV, 

where the “ 0 ” subscript indicates that the quantity is evaluated on the background scalar 
field configuration; e.g. := We further specialize the model by choosing the 

matrix 7 ,^^ to be independent of the background field 

7 ^'® := = constant, 

but retain the freedom that may contain off-diagonal elements. The action then reads 

^2 = I (v^/aV^/s + KIaIb) dV + 5^pp, (4) 

where 

dSp^ = -j q^fAdV. 

By performing an integration by parts we see that the action can be written in standard 
quadratic form by introducing an operator given by 

(□ - V ’’), 

such that 

^2 = -^ y fAO^^fB dV + <5^pp. (5) 

The inverse of the operator defines the propagator of the theory. The propagator, along 
with higher n-point functions, are obtained via the moments of the generating functional 

where is an arbitrary external source. The smallness of the particle mo ~ iJa ^ 1 permits 
a Taylor series expansion for e*'^pp, after which ^ becomes 

nj] = j (^l + t6S^^ + '^5Sl^ + ---^e^d^^fAdv_ ( 7 ) 

Replacing /a by the variational derivative —i5/5j^ taken at = 0, we reduce the above 
path integral to a Gaussian in curved spacetime. Although in a general curved background 
computing the Gaussian path integral is not tractable, it is manageable in fiat spacetime. 
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and we may exploit the equivalence principle to design a perturbation series based on local 
flatness about a point. To do so, we work in a local convex neighbourhood ^{xq) centred 
at a fixed point Xq, and expand deviations from flat spacetime in powers of the coordinate 
distance from xq. In this paper we adopt Riemann normal coordinates (RNC) which we 
denote by x°‘. The Riemann normal coordinates of a point x with respect to the fixed base 
point Xq and tetrad e°‘^{xo) are defined by the relation 

x°‘ := -e°‘^{xo)a^^{x,Xo), 

where we have introduced Synge’s bi-vector a'^{x,xo) which runs tangent to the unique 
geodesic connecting xq and a point x in ^{xq). The RNC expansion is valid provided that 
the metric varies slowly on scales of order the size of the convex normal neighbourhood. 
A full discussion of bitensors including Synge’s function can be found in the review article 
[21]. As the held operator is built solely from background quantities, it admits the 
expansion 

O'^^ix) = O^^(a;o) -h -h O'^jl{xo)x^x'' H-, 

:= ( 8 ) 

where denotes the flat spacetime kinetic operator given by 

= ■y-’® /a /b, (9) 

which yields the principal part of the held equations upon variation. The quantity 
incorporates the curvature corrections to the bulk propagator in the local neighbourhood of 
Xq. For general fields, will also contain contributions from 0^^{xo) due the presence 
of background curvature coupling directly to the perturbed fields, as is the case of vacuum 
gravity through terms such as R^'^^^haph^s- Such terms persist independently of any ex¬ 
pansions the background metric might have. The metric determinant found in the volume 
element dV = yf—gd'^x also contributes curvature terms in the RNC expansion and these 
must treated with care. The metric determinant in the volume element of exp (— J dVj^f^) 
can be eliminated by shifting the external current This shift is of no con¬ 

sequence because we set the current to zero after variation. To take care of the metric 
determinant in the bulk action coupling to 0'^^{xq), we use the fact that g is expanded 
around its flat spacetime unit value and we absorb the higher-order terms into a redefinition 



of As is represented by a series in powers of its contribution to the bulk 

phase can be written as a series analogous to that involving the internal source in 

Eq. (7) with again being replaced by a variational derivative with respect to the current 
j^. Explicitly, we hnd that to second-order in the derivative expansion, the quantity S’ab 
reads 


= fA fA- (10) 

After performing the Gaussian integration 


exp 


'2 / fA^^^fBdv + i / j^fAdv 


= =:^b1, (11) 


7-B' 


where dv is the volume element of Minkowski spacetime and a ■ b := J a{x)h{y) dvxdvy, we 
hnd 


= 1 - yy DABiz,z')dTdT'+ jj drdr' j DAc{z,x)S’^'^{x)DDBix,z')dv-\ - 

( 12 ) 

after functional differentiation. We then Legendre transform the generating functional with 
respect to the world line coordinates to obtain a formal expression for the effective action 

rW = Rpp - J dr J dT'GAB{z,z'), (13) 

where we have used the notation Gab for the complete curved spacetime Green function; 


i.e., 


Gab — Dab — i j dv Dac^^^D^b + 0{<S’^). 

Varying with respect to the world line coordinate gives the equation of motion 
m(r)a“ = q^w'^Np^A + p f dr'G^^z, z!), 

J —oo 

where 


(14) 


(15) 


projects on the spatial plane orthogonal to the world line. The equation of motion must be 
expressed in terms of the retarded Green function to satisfy causality. Although at this stage 
the retarded propagator does not appear to arise naturally from the effective action, which 
is time-symmetric under interchange of r and r', we can ensure that the Green function 
is retarded by adopting the causal action formalism to break the time symmetry resulting 
from a naive application of Hamilton’s Principle [32]. The causal action formulation and 
the conditions required to obtain the retarded propagator will be discussed in a later section 
(Sec IV). 
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B. Quasi-local expansion 


As it is written above the equation of motion still suffers from a divergence on the world 
line due to the singularity present in at the simultaneous point z' = 2 :. In order to 
proceed we must regularize the expression, which will allow us to separate the equation of 
motion into a sum of hnite local terms and a tail integral over the particle’s past. The 
method of choice for regularization is dimensional regularization, which is well adapted to 
the EFT formulation. Dimensional regularization requires that we work with a momentum 
space representation for the flat spacetime retarded propagator To regularize the 

equation of motion, we must hrst determine its local behaviour. Recall that we have a local 
expansion in RNC for the bulk curvature operator around a base point Xq as given to 
second order in Eq. (10). Since xq is arbitrary we can choose it to coincide with the 

coordinate of the world line at the “current time”, so that the origin of the RNC system is 
comoving with the particle. Likewise we choose x" to be the RNC of What’s left is to 

translate the Riemann coordinates, which are expressed as a series in the geodesic distance 
between the two points on the world line through Synge’s bivector (Tq, into the quantities 
dehned on the accelerated world line such as the four velocity and its derivatives. Following 
[36], we translate between quantities written in terms of the geodesic and the accelerated 
word line connecting the two points x and Xq by introducing a small parameter t dehned as 
the proper time difference between the current time r and some time in the very near past 
t'. Once the quasi-local expansion in terms of t is hnished, the expression for the force in 
terms of a tail integral and an integral over Fourier modes can be manipulated into a form 
which lends itself to dimensional regularization resulting in a hnite expression for the force 
including local terms. 

We now compute the local force 

j dr'D%{z,z'), (16) 

due to the leading-order term in the RNC expansion of the Green function found in Eq. 
(14), which diagrammatically corresponds to computing graphs without insertions of 
In Eq. (16), is evaluated at the present time r and the retarded propagator is expressed 
in RNC centred at z{t) as the Fourier integral 

D%(0,x)=Aab [ (17) 

Jk,ret ^ 
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where Kab is defined by the relation Kac^~^— ^ab, is the Minkowski space inner- 
product of wave-vectors -|- /c^, and the subscripted integral corresponds 

to a d-dimensional Fourier integral using the retarded contour in the ca-plane. The u integral 
has two poles sitting at plus and minus the magnitude of the spatial vector k, which forces 
us to choose a contour for computing the residues. The subscript ret indicates that we have 
chosen a contour that passes over the poles and closes at Imca —)■ —oo corresponding to the 
condition r > r'. Then, using the definition of RNC we make the substitution x'^{z') = 
—e'^p{z)a^{z,z'), and orient the tetrad such that e^Q = This allows the exponential in 
k^AB to t)e expressed in terms of aa{z, z') as defined on the geodesic connecting and z' as 
gifc o- terms of the accelerated world line, the geodesic bivector has the expansion 

ao,{z, z') = -tUa - H-, (18) 


where Oq, is the covariant acceleration Oq, = iia ■= Dua/dr. Using the above expansion 
for Synge’s bivector, along with the definition of RNC, we find that the exponential in the 
Green function takes the expanded form 


^ik’X _ 


ika'^ 


fn+l jjr 

^ (n -I- 1)! dr’ 


-I- 


(19) 


where the numbered/latin indices indicate frame components; e.g., /c“ = /c^Cq-|-/ c*e". Chang¬ 
ing integration variables from r' to t using t = r' — r, we find that all integrations necessary 
for the local expansion of the force involve evaluating a single master integral 




dt 




e{kyh,---h^, 


'k Jk^ {k^ — {k^yA 

■2n-p (d + p + g-2 m-l)! 


r (^) r(m) 


(g + 1)!! 


dn^d+q+p—2m—l 


( 20 ) 


where is a sum over the product of Kronecker deltas with all possible permutations of 
q indices, + perms. The above integral is an extension of the master 

integral used in [38], which was previously specialized to the single m = 2 case. The integral 
is straightforwardly performed after making a few observations about its form. The first is 
that odd powers of spatial momenta yield zero by rotational invariance. Secondly, powers of 
t can be written as derivatives with respect to k^ after integrating by parts. After integration 
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by parts, the t integral becomes a delta function 5{k^), and the spatial k integral is reduced 
to a standard form that can be performed via the Beta function 

/■“ fc"-* _ .,,_,„ r(V2)r(m-V2) 

Ja (k^ + v^)"' 2r(m) 

The result of the Euclidean integral, along with the distribution S{k^) coming from the t 
integral, results in the Kronecker delta 6 n,d+q+p- 2 m-i whose arguments are determined by 
the exponent of k^ which must be zero otherwise the whole expression is zero [38]. The 
upshot of this is that t power counts like l/|/c| and k^ doesn’t need to be power counted 
when determining the high energy behaviour of the integral. 

Using the above expansions in conjunction with the master integral we hnd just a single 
non-zero term in the series (19) involving the integral Jij(l,3,0,2) = ^Sij leading to the 
local self-force expression 

fa = Waisa^. 

4:71 6 

Notice that in the limit A = B = 1, the local self-force above reduces to the expected result 
for a single charge [29]. 

Let us now examine the effect of curvature terms in the bulk action on the local self¬ 
force. The presence of bulk curvature entering through the quantity S’ab causes a focusing 
of null rays and the detailed history dependent force which results is generally so complicated 
that it must be evaluated numerically or left as a formal expression. The local nature of 
the force in curved spacetime, however, can be established in the normal neighbourhood of 
the particle where unique geodesics exist. The starting point for the computation is the 
RNC representation of given in Eq. (10). Power counting, we hnd that the high energy 
behaviour of the potential term in four dimensions is 0(A“^) for a cutoff A ~ (size of body)“^ 
and will not lead to a local term. To second-order in the RNC expansion we hnd that the 
propagator (14) reads 


Gab{Xo,x) = J y~'ll‘SADAc(X(„ y)S:,DoB{ll,x) + --- 


( 21 ) 

where the Riemann tensor is evaluated at the base point Xq at which the Green function 

Dab has the form (17). We hnd that the term f j^y'^y^daDAcdiiDDBiy^x') vanishes 

algebraically and the only contributing part at second-order in background derivatives reads 

1 /■ 17 

j y^y^dpDAcd^^DoB = --Aab {Rk‘^ - 2R^yk^k^) , (22) 
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where we have eliminated using the relation After inserting the 

expression above into Eq (15), we find the local self-force containing curvature terms is 
given by | WayR^'^u^ as expected from the known result of the self-force on a 

single charge. 


C. Equation of motion 

Assembling the local expressions, we hnd that the explicit self-force equation of motion 
reads 


^ ^ / dT'G%{z, z'). (23) 


III. MULTIFIELD FORMALISM 


We now consider a more general scenario in which a small non-rotating but otherwise 
arbitrary object sources a set of perturbations of a background non-vacuum spacetime. 
The background spacetime, ^g, is dehned by a smooth metric satisfying the non-vacuum 
Einstein held equations in the absence of the small body (or a suitable limit where the 
mass, charge, and size of the body go to zero uniformly). We characterize the small body 
by a set of internal sources, which we collectively denote by a set of local quantities 
characterizing the multipoles of the body. In the multi-scalar held toy model, the source 
took the form = —g^ f drS^x, z). The source generates a perturbation of the background 
helds, which we denote collectively by the held multiplet ip a- Eventually, we will consider 
the set of sources = {p, where p is the scalar charge density, is the electric 

current, and is the stress energy of the world line, along with the corresponding set 
of held perturbations ip a = hap}- With this notation, the discussion very closely 

parallels the treatment of the multi-scalar held toy model. 

We begin by constructing a bulk action quadratic in the helds which we add to the 
phenomenological action characterizing the extended body as a point particle with a set of 
multipole couplings. We construct the quadratic action by applying a perturbation operator 
to the full non-linear Lagrangian of the theory and integrating the result over the background 
manifold. For our purposes we dehne the n-th order perturbation of a quantity Q as 5^ Q ■.= 
In a local coordinate chart of the background manifold with orientation 60123 = 1 
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having four-volume weight \/—g, the second-order perturbation of the action reads 


^2 


1 

2 


dv - 


fi^tl^AdV =: + 5^pp, 


( 24 ) 


where dV is the invariant volume of the coordinate chart. As we did with the toy-model, we 
separate into a part which is invertible in flat spacetime, and a curvature oper¬ 

ator which we treat perturbatively in the local neighbourhood of the body. In writing 
Eq. (24) we’ve assumed that the held variations (and their hrst derivatives) vanish sufh- 
ciently fast at large distances to ensure boundary terms can be discarded when performing 
integration by parts. 

As we did with the multi-scalar toy model, we start constructing the effective action by 
introducing a generating functional which produces all the possible held correlations of the 
theory. To build the generating functional we introduce external sources which couple 
linearly to the helds. The generating functional as a functional of the external current 
coupled to the held reads 

^[J] = J (25) 

As fi^ is perturbatively small, we express as a series expansion, after which ^ becomes 


S’lJ]= j (l-i j dVdV + ■ ■ A e‘d , 


(26) 

just as for the toy-model. The right-hand side of Eq. (26) is simplihed by replacing iIja by 
the variational derivative —iS/SJ^ evaluated at = 0. With all powers of ip removed we 
are left with a Gaussian integral in curved spacetime. In the toy model we simply evaluated 
the Gaussian and proceeded from there. However, performing the Gaussian integral requires 
to be invertible. If tp contains helds with gauge symmetries, the inversion is problematic 
because becomes a singular matrix and the path integral diverges due to over counting 
physically equivalent held conhgurations. To remedy this, we adopt the Faddeev-Popov 
prescription which involves adding an additional quadratic contribution to S 2 given by 


= - j GAG^dV (27) 

where G^ is a gauge-hxing constraint acting on the helds. We state without proof that the 
Faddeev-Popov prescription gives a satisfactory result without introducing spurious gauge 
modes. Further details regarding the Fadeev-Popov prescription will be elaborated in the 
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appendix. We absorb the extra terms resulting from the gauge fixing procedure, which 
involves several integrations by parts, into a redehnition of , which we will call 
The gauge-fixed operator is invertible (at least in a perturbative sense). 

Following the steps in the toy model, we again work in the normal neighbourhood of a 
point and write the inverse Green function as a coordinate expansion where the leading 
order term is the invertible in flat spacetime and curvature corrections are contained 
in the quantity With this, the curved spacetime propagator takes the form of (14) as 

it enters the generating functional 


^[J] = 1 - 


1 


dVdV'] Zo[J]. 


( 28 ) 


So far we have left the internal source arbitrary, merely requiring that it be small 
in magnitude relative to the scale set by the background. To proceed further we consider 
the specihc case for which is an effective point-particle world line. For this description 
to make sense in the context of effective held theory, we should also add additional higher- 
order couplings involving non-minimal “operators” which encode information about the 
internal structure of the body. However, it has been shown that these higher order terms 
are irrelevant for a hrst order self-force calculation, and so we ignore such terms here [35]. 
A point source amounts to the choice of currents 



where specihes the type of current in question. Since we have now introduced a point 
particle degree of freedom into the theory, we must also introduce an external source 
coupling to it, and from it we build a generating functional for the entire system 

Z[J,C] = J (29) 

where the subscript pp serves to indicate that the point particle source is to be used for 
/i"^. Although we have explicitly introduced a coordinate dependence into the generating 
functional, this coordinate dependence will ultimately cancel in the effective action. The 
connected portion of the generating functional is obtained via a logarithm. Performing a 
Legendre transform of the connected part with respect to the external current coupled to 
the world line yields the effective action 

T[z\ = -i\nZ - [ Cz^dV, (30) 
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which generates all one-particle irreducible connected diagrams. Quantum fluctuations of 
the world line are so highly suppressed for macroscopic astrophysical bodies that we may 
treat z a.s a classical word line. The upshot of this is that the integration over all world line 
configurations is easily performed using the saddle-point approximation. With the addition 
of the zeroth order point-particle action which is responsible for accelerated motion in 
the background, we arrive at the effective action for our problem 

dTctT J dv D^D£)B '-!-■■■ 

:=‘S'pp - ^ jj drdr' s^Gab'S^', 

where is evaluated at a bulk point in Riemann normal coordinates. Although the series 
for Gab proliferates indefinitely, on the basis of dimensional analysis only the first two terms 
which we have explicitly written down contribute a local self-force. All higher order terms 
(0(^^) and above) are contributions to the tail which we will eventually represent in terms 
of a formal expression. 

IV. CAUSAL EQUATIONS OF MOTION 

To derive an expression for the self-force using an action principle, one must adopt a 
formalism where the effective action describes causal evolution under the effects of radiation 
reaction. In the quantum held theory of non-equilibrium systems [39] this formalism is 
commonly referred to as the IN-IN construction or closed-time-path (CTP) formalism. It 
differs from the more common IN-OUT formalism, where scattering boundary conditions 
are placed at the remote future and remote past and the propagator of excitations is the 
Feynman propagator which is symmetric on interchange of spacetime event points. In the 
IN-IN formalism the IN states are prescribed and the future states are determined by the 
evolution of the IN states under the action of the generator of time translations. A classical 
formulation of this has been introduced by Galley under the moniker Causal Action Principle 
(CAP) [32, 33]. Under the CAP construction one breaks the time reversal symmetry by 
introducing two copies of each configuration variable in the action. One then imposes the 
conditions that both variables are stationary at some initial time (zero variation) and that 
at some final time, which may be in the remote future, the two configurations are equal but 
their values are left unspecified. Causal equations of motion are obtained by varying the 



PM = 


pp 


- JJ drdr s"^Dab'S^' + 
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paths and then setting them equal to each other. 

The CAP action for this theory is 

S = Si - S 2 := 5'bulk[^A,l] - *S'bulk[^A,2] + Spp[z^,ijA,i] - Spp[z2,iJA,2], (32) 

where Si is the action obtained by evaluating all quantities on the “1” conhguration and 
S 2 is the action obtained by evaluating all quantities on the “2” configuration. The 1-2 
index structure of the RHS above suggests the introduction of a 2-d space with CAP indices 
a,b,c,.. = 1, 2 with Lorentzian metric 

C'a, = C''“ = diag(l,-l). (33) 


To produce the generating functional we couple the field configurations and world lines to 
external currents. With these considerations, the generating functional for IN-IN correlation 
functions in this theory takes the form 


z[JX] 


Qiz Slip exp 


'CTP 


iS + iC 


ab 


Ca^ba + 



(34) 


The functional measure for the path integral is dehned by the conditions that the configu¬ 
rations are equal to some unspecihed value at a time t = tf in the remote future and that 
their values are specihed on a spatial slice S(tj) in the remote past by suitable initial data 


= 


'CTP 


/E(q) 






X 




^Ipl 






S'lp-i, 


(35) 


As we are only interested in classical trajectories of the particle, the measure for the particle’s 
paths can be ignored. The propagators for the theory are given by functional differentiation 
with respect to the external currents 


D%,{x,x') 


1 PZ 


(36) 


Notice that the two-point function is now a matrix-valued quantity, each element of which 
corresponds to propagation on a different branch of the CTP. The matrix for a scalar 
held has the form 


T>“'’(a;o, x) 


Df{xo,x) -D_{xo,x) 
~Df{xo,x) Dd{xo,x) 


(37) 
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where Dp, -Dd, and D_/^ are the Feynman, Dyson, and negative/positive Wightman two- 
point functions, respectively [35]. Using the saddle-point method, the effective action is 
readily found to be 


hcApi^j = // dTadrl s^D%,sf+^- // dTadrl / dv D%,sf + 


(38) 

where all CAP indices are summed over. Before taking the variation to obtain the equations 
of motion for the particle, it is useful to introduce average and difference coordinates 


:= 


(39) 


The benefit of the above coordinates is that, after variation when we impose the equality 
condition Zi=Z 2 -i z°!_ vanishes and 2 :" can immediately be identihed as the physical path. 
The equations of motion governing the physical path z^ are found by varying TcApi^j with 
respect to 2 ;“ by demanding that the physical path is a stationary point 5Vqxp[z\I dz- = 0. 
The causal equation of motion is found to be 


0 = 




z _=0 


+ 

2 ' 


dT'D%,s^ --F/ 


dv dr' + ■■■ , (40) 


where the tensor F^{t) decomposes into gradient (acting on the Green function) and accel¬ 
eration components 

n ■-= 

and has the useful property Dab'F^u'^ = 0. As all world line quantities in the above 
equation are evaluated on the averaged (physical) conhguration, we shall omit the “-I-” 
index from now on. 


V. QUASI-LOCAL EXPANSION 

Although we have succeeded in obtaining a causal equation of motion by extremizing the 
effective action, the expression is ill-dehned because it contains a divergence at the position 
of the particle coming from the singular nature of the retarded propagator D'’®* ~ (5 ((t) on 
the past lightcone. In this section we present a local expansion around the accelerated world 
line for an arbitrary source that extends the local treatment of the multi-scalar field to 
helds with scalar, vector, or tensor representations. 
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A. Local multi-field self-force I: ignoring curvature contributions from 

We first compute the terms in the equation of motion which do not result from curvature 
contributions coming from insertions of . The procedure we apply here generalizes to 
the terms with curvature as well, but the discussion is simplihed without their presence. 
The force which persists in the absence of spacetime curvature under the definition 

Ut) j dr', ( 41 ) 

where is evaluated at the present time r and the retarded propagator is represented in 
RNC with origin at z{t) as a Fourier integral 

, (42) 

Jk,ret ^ 

and where Aab is dehned by the relation = ^ab- An example of the type 

of force produced by fQ(T) is the Abraham-Lorentz-Dirac force describing the radiation 
damping effect on charged particles in flat spacetime. 

The local expansion of Eq. (41) proceeds in the same way as in the toy-model; we write 
the Green function D^b and source as a series in powers of Synge’s bivector on the 
geodesic connecting 2 ; and z' and then relate these quantities to the quantities dehned on 
the accelerated world line using the expansion of in Eq. (18). This leads to an expression 
with several terms all involving the integral (20), which we now examine. Inspecting (20) 
with the hat spacetime propagator D{k) ~ /c which implies m = 1, we see that the 
integral yields zero in d = 4 unless the condition n = g + p + lis satished. 

We begin by separating the integral expressions according to whether they involve either 
or constituting F^. First, we consider the contribution. Investigat¬ 

ing several terms in the local expansion, we hnd that the condition n = g + p + lis satished 
by one term only. The reason is that after removing pieces containing odd powers of spatial 
vectors k\ which average to zero, the leading terms are (1 — i(t^/3\) k^ao){s^ + ts^ + ••■)) 
and since g = 0, the only term remaining is ts^. From this we get a contribution of the 
form A\i^a^s^. Now, as must yield a reparameterization invariant action, it must be 
constructed from the four-velocity alone. It follows that A^^a^s^ is quadratic in the accel¬ 
eration. Higher order terms in /c“, including those with q > 2 which involve products of sums 
containing multiple powers of t for a single increase in the power of k^, cannot satisfy n = 
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p+q+ 1 and yield zero. We now move to the part of F^. The presence of the covari¬ 

ant derivative introduces an additional positive power of /c". Power counting the exponents 

in the expansion {k^u^ + — /2\) kja^ —i{t^/3\) (k^do + Fdi) -\ - ){s^ + ts^ H -), 

we see that three terms can satisfy the non-zero condition n = q + p + 1 with q even: 
two proportional to a derivative of the acceleration and the other one proportional to the 
acceleration squared. 

Collecting the above results, we find that the self-force in flat spacetime with a general 
source is 








1 D 

2 dr 


s^u^) 


+ 




(43) 


The equation above can be evaluated immediately once is determined, and so it provides 
a very efficient route from the variation (40) to a hnite part of the self-force in non-vacuum 
theories. We now move to deriving a similar expression that includes curvature corrections. 


B. Local multi-field self-force II: curvature coutributious iuvolviug 

The focus of this section is the local self-force due to the curvature expansion of the 
equation of motion (40) in powers of insertions. First, we introduce the quantity 

fc - j [ ( 44 ) 

which dehnes the leading-order, curvature-influenced local force. At this order, the curvature 
operator only contains single powers of the curvature tensors and they appear undifferenti¬ 
ated. 

To proceed, we require a slight extension of the procedure used in the previous section 
to accommodate the additional coordinate expansion of the bulk curvature operator = 
where y" is the RNC of a bulk point in the normal neighbourhood of the origin. 
The introduction of uppercase Latin letters with carets facilitates the use of multi-indices 
L = Pi ■ ■ ■ /i£. In order to use the master integral expression (20) we must recast fa into the 
appropriate form. To do this, we need to perform the volume integral. The volume integral 
in RNC centred at z{t) has the form J where, as before, x°‘ is the RNC 

of the “past” point Rewriting as and doing ^ integrations by parts 
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reduces the volume integral to a delta function S{p — k), leaving a single Fourier integral. 
Transforming back to global world line coordinates we obtain 


fa = — ^F^^ac^db J J 


La 




dki k"^ 


+ F 




^CD 
1 


F' 


V dki F 


dki F 


k^ 

d^ 


(45) 


yCD _ 

dkpF 


^ik-asB' ^ 


where and represent coefficients of terms in the bulk curvature operator 

one, and zero derivatives, respectively. We therefore refer to these coefficients 
as the kinetic, angular momentum, and mass tensors. The tensor has the symmetry 

On dimensional gronnds, begin at £ = 2, £ = 1, 

£ = 0, respectively. For reference, in the case of the multi-scalar toy-held we had = 


{\RaMS - \FsRap) and 

Again, we analyze the self-force expression in terms of the acceleration and gradient 

R^af^^ pieces that constitnte the tensor The part involving gives zero due to 

the additional powers oil/k coming from products of the hat spacetime propagators. Its 
leading order piece, which contains the integral 1(2, 0, 0, 0) is convergent for large k like A“^. 
All sub-leading terms in the sum contribute additional k~^ powers and serve to make the 
result more convergent. The gradient component does contribute local terms to the 

equation of motion, however. The hrst of these terms is derived from the RNC expansion of 
the metric determinant which enters as the £ = 2 component of the tensor coefficient 
Using the power connting based on onr master integral expression, we see that mnltipoles 
higher than £ = 2 converge rapidly in the local expansion. For £ = 2, we obtain the integral 

^ j dtj^^ (^2r|F^^F'^ + pf^’^F - k^F^°\ 

which involves combinations of the integrals /jj(3,0,1,2) = —i/327iSij, /ijfci(4,0,1,4) = 
i/192nSijki, and 1(4,0,3,2) = —i/192nSij. Althongh there are a large nnmber of per¬ 
mutations involved, the expression is simplified somewhat by the orthogonality condition 
= 0 and the need for an even nnmber of spatial k vectors. For £ = 1, we hnd the 


integral 



which has non-zero contributions conpling to the vector coefficient from the integrals 


/jj(3,0,l,2) = —i/32Ti5ij and 1(2, 0,1,0) = —i/Sn. Lastly, we have the £ = 0 expression 
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which is simply the integral /(2,0,1,0) multiplied by a single four-velocity vector. All the 
above integrals may be performed independently of the representation of the held, but the 
hnal expressions are rather lengthy and considerable simplihcations occur when a specihc 
choice is made for the theory-dependent quantities and 5"^^. Considering their length, 
we delay their appearance until we specify the held representations in our discussion of the 
scalarvac and electrovac spacetimes. 

We end this section with a formal expression 


/<.(t) = 


■'^B-Vh,nh 


ACJ-^DB 


n^CD j I c/>CD T y/CD j 

l3 ' “T -'/S 


R 


( 46 ) 


for the local self-force due to curvature terms in the bulk action as a local function of 
proper time. Although Eq.(46) appears more complicated than its hat spacetime counterpart 
Eq. (43), as we shall see later when we consider the scalarvac and electrovac space times, 
the quantities and are easily obtained from the second variation of the 

bulk action. Additionally, most of the contractions involving the Riemann tensor with the 
integrals and J/j give zero, leading to roughly the same number of terms as (43). 


VI. SCALARVAC SELF-FORCE 

A. Perturbed scalarvac action 

We now move to our hrst concrete example: the perturbed scalarvac spacetime. Here 
we consider a background spacetime with metric gap and a minimally coupled background 
scalar held <h with action given by 

^buik = V„4>V^4> + 2V(4>)) dV. (47) 

In four spacetime dimensions, the background helds evolve dynamically according to the 
held equations 

Rap = Tap — \gapT^ (48) 

and 

□$ _ W($) = 0 (49) 
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where = dV/d^. The stress-energy tensor of the background is given by 


Tal3 = 


2 6S, 


matter 




( 50 ) 


which implies = V+ g^^V (*h) for the Ricci tensor. The background spacetime is 
perturbed by the presence of a small body, which at signihcant distances from its centre looks 
like a point particle with mass m, scalar charge q, and higher order multipolar couplings 
which we will ignore here. The particle moves on a world line 7 described by parametric 
relations z^{t), where r is proper time in the background spacetime. We denote the four- 
velocity in the background spacetime as m" and normalize it in the background metric such 
that gapu'^u^ = —1. We add to our background action the following point particle interaction 
term 


Rpp = -m 


dr + q J (ir<h(r). (51) 

The disturbance introduced by the particle sources a perturbation of the helds dehned 
as the difference between the full non-linear helds, if and and their background values 


/ = (p - $ 

ha/B ga)S- 


The two perturbed helds are packaged using our multi-index notation into the doublet 
BjjA ■= {/, ha/s}- The variations ddr = —^h^uU^u'^dr and 6^ = f couple to a source doublet 
drS^i^x, z)s^{z) given by 

:= {q, l/2(m — q^)u^u'^}. 


We shall also require the hrst two derivatives of the source 

= { 0 , {m — g<h) 

and 

= {0, (m — g<h) -I- a^a^) — ^q^u^u" — 

in the expansion of the world line. 

To apply the multiheld formalism in Sec. Ill, we expand the action to second-order in 
the held perturbations and introduce a gauge. This allows us to compute the propagators 
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associated with the helds 'ipA- After removing terms proportional to the background field 
equations, we hnd that the second-order variation of the bulk action reads 

5'^buik = j (^-|V«/V“/ - h,s (52) 

+ dV, 

where 

^a/37-5 = 1 + g^^g^^ - g^^g^^) 

is the identity on the space of symmetric, trace-free, second-rank tensors. Gauge hxing is 
implemented via the Faddeev-Popov prescription. For this we choose the gauge condition 
Ga = I {La - 2/V„<h), where La = [hag - \gaph) is the Lorenz gauge vector. The 
beneht of choosing this gauge is that it exchanges the derivative coupling between the scalar 
held and the metric perturbation for an algebraic coupling. After simplihcation we hnd 

i"Sb„ik = -\ j [f (-a+ M)f + \h„0 / )dV, 

(53) 

where 

^ := P"($) + 2V„$V“$, 

j^fMiyag 2 _ 2 P($), 

^a/3 _ 2V“V^<h. 

After second variational diherentiation of the bulk action we hnd the gauge-hxed held op¬ 
erator, which has the matrix form 

-D + 

Notice that in hat spacetime the matrix has the diagonal form 

o,'‘® = diag(^-n/,-lAA‘n,). 

where □/ is the hat spacetime wave operator g^''dfj,di, and has the index / to indicate 

that the hat spacetime metric is to be used in the expression for To satisfy the 

equation O^^Dbc = spacetime Green function must equal 

Dab = diag ^ , 
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where PajS'^s is defined such that | [S^Sg + S^Sgj. In d spacetime dimensions 

pjp^s has the form = 1 _ Note that in four spacetime 

dimensions PajSyS and are equal. Additionally, we have the causality requirement 

that 0“^ be chosen according to retarded boundary conditions, as dictated by the CAP 
construction. 


B. 


Local scalarvac self-force I; expansion for vanishing 


Here we set out to compute the quasi-local expansion of the scalarvac self-force in fiat 
spacetime according to Eq. (43). For this we need to determine the quantity which 
results from varying the effective action 




CAP 


drdr'= 2 j dr' 






5s^ 


rgret 

^AB' 




(55) 


The scalar source s = q, which has vanishing time derivative, leaves the simple expression 
f dr' (tta + D[z, z'), from which we readily read off Fa = 2q{aa + ^). Thus, 

Aap = 2 qgai3 and Bap = 2 qwa(3- The tensor element of F^ is complicated by the presence 
of the time dependent source (m — q^)uaUi 3 . To simplify matters slightly it is useful to 
introduce a time-dependent mass given by 


m := m — g<h. 


After the variation 6 s^/Su^, we find the expression ^ — \u^u''aa) m{T)D^yp 

where the derivative acts on both the Green function and the mass function, and 


10^'" = -u^u 'w" 
2 




u ’'u . 


Using the background equation of motion moa = qw^Vp^, we write pm as 

— ^mu^u'^aa + qwa^u^^^ and lump it into the piece of that doesn’t include 

a derivative operator acting on the Green function. We proceed by separating FJ^^ into 
pieces which have a derivative on the background scalar field instead of the original sepa¬ 
ration involving the acceleration, and those which act on the Green function. With this. 
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we have = —q coupling to a derivative acting on the back¬ 
ground scalar held, and coupling to a derivative on Dab- Substituting 

these values and the expressions for the derivatives of the sources into Eq. (43), we arrive 
at the local contribution to the self-force with vanishing 

47rf„ = + Tgrua^®, (56) 

with the overall factor of 47r deviation from the result of [1] being due to normalization differ¬ 
ences in the action. Notice that we recover the well-known gravitational anti-damping and 
scalar damping terms, —and respectively. We also hnd a new term, 7gmaQ,<i>, 

due to the additional coupling between helds. To reduce the order of this equation we use 
the background equation of motion and its derivative in the orthogonal 

direction ho = . After simplihcation we hnd 

4irf„ = (A _ 1 j + img (-11 + A'j (57) 

for the order-reduced local (hat) self-force in terms of derivatives of the background scalar 
held. 

C. Local scalarvac self-force II: contributions from 

Here, we calculate the local self-force resulting from curvature terms in the bulk for the 
scalarvac theory. These terms are encoded in the curved spacetime propagator through 
the tensor In a general curved spacetime, the global propagator is not exactly solv¬ 

able but it can be represented as a perturbative expansion in a small region of space- 
time known as the convex normal neighbourhood. One approach, due to Hadamard, is 
to assume a distributional ansatz with the appropriate causal structure Gab'{x,x') = 
UAB'{.x,x')5{a) + VAB'{x,x')9{—a), for which the smooth biscalar coefficients U and V are 
determined recursively using the held equations. Using the Hadamard form, one can apply 
the Detweiler-Whiting regularization method [30] to obtain the helds responsible for the 
self-force. This was the approach used in [1]. Here we do not employ the Hadamard form, 
but instead we build the curved spacetime propagator Gab as a series involving the hat 
spacetime propagator Dab and corrections coming from the curvature terms ^ab in the 
bulk action just as we did with the toy-model (21). 
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To begin, we break up the bulk action into three parts 


<^^>S'buik — Sh + Sf + s„ 


consisting of 

Sf = -\f dv{g^^djd^f + (l""(<h) + /'), 

Srm. = ~J dVKp {2d^d^ + f, (58) 



where we have used the coordinate dehnition of the covariant derivative Vfj,ha /3 = d^hag — 
T^afidup — In writing the coordinate form above, we have anticipated the use of 

Riemann normal coordinates which we now reintroduce. In the local convex neighborhood 
of xq, the scale TZ at which the metric changes is much larger than the typical length of 
the coordinate to a bulk point y". RNC therefore facilitate an expansion of background 
quantities around their local values as a power series in the quantity y°‘da ~ r/7^, where r 
is the magnitude of the geodesic length connecting the two points Xq and y. We will refer to 
this expansion as the near-point expansion and terms with n derivatives of the background 
helds will enter at n-th order in the expansion. Using the notation of Sec. Ill, we introduce 
a bulk interaction Lagrangian in the local convex normal neighbourhood given by 


Sint — 


dv, 




1=2 




c>AB,a L—1 


AB „.L-2„ 


(59) 


where the coefficients and are evaluated at the base point xq. 

For the hrst-order self-force, we only need to compute the next-to-leading order curvature 
corrections to the propagator, and so we drop terms in the near-point expansion of S’ab with 
more than two derivatives. Reading off the quantities from equations Eq. (54) and Eq. (58), 
and dropping terms third order and higher we hnd 

0 

0 
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0 


0 


0 


vX 


and 


2do,dp^ 

Using the RNC expansion of the backgronnd metric 

9o/j = vcfi-^R<,-,gsv''v‘ + ■■■, g°^ =+ j-R“U y''y‘ + ■ ■ 


(60) 


and related quantities 

9 = ^- ^Ra0y'^y^ + • • •, + • • •, 

we hnd that the kinetic, angular momentum, and mass vertices in RNC read 

+ \r“;„ 


1 


+ 5-)'“' {RW’^ + Ry',v‘^‘ - + nto”'’)). 


= 0 , 

o 

^a/3'i5,fj, _ j^a^ri-yj^S) 

3 

^ = u" + 29„d>(9"$, 

^a/3 _ 4di»dy)<i>^ 

— ]Yfial3)iy5) ^ 


(uT)y^ 


(61) 


(62a) 


(62b) 

(62c) 

(62d) 

(62e) 

(62f) 

(62g) 


where we have introduced the shorthand notation Ay = A^y^ to indicate contraction with 
bulk Riemann coordinates. We now replace powers of the Riemann coordinate with deriva¬ 
tives with respect to /c“ to put the equation in the form of Eq. (45). 

Evaluating Eq. (46), we hnd that the part of containing the term 

R'^a dk^^dkp W vanishes by virtue of the symmetries of the Riemann tensor. We also hnd 

that the entire integral is canceled tensorially by the second line of 
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in (62b) when coupled to . Inserting the remaining expressions for the bulk vertices 

Eq. (62) into Eq. (46), eliminating terms by orthogonality, and employing the integral con¬ 
dition n = p + q + d — 2m — 1, we hnd non-zero results from the metric determinants in 
and mixing tensor and the Ricci curvature constituents of In 

the formulation using the Hadamard form, the part of the form corresponding to the metric 
determinant pieces of and would be the van Vleck determinant in the direct 

piece of the curved spacetime propagator. 

After performing the integrations we find that the scalar held kinetic term containing the 
metric determinant gives the local term the tensor kinetic term and the tensor 

mass term give the contributions '^Rapu^ and —2m^V, respectively. The term in 
the local force due to the mixing is 2mgVoVy3<h . Combining everything, we hnd that 
the local self-force resulting from bulk curvature insertions reads 

An fa = (f‘)Ri 3 'yU^ — 2m^u^ V/ 3 <hV.y<h -h 2gm . (63) 


D. Scalarvac equation of motion 


Amalgamating the local results just derived with the tail terms coming from the variation 
of (55) with the full non-diagonal Green function Gab' we arrive at the covariant equation 


1 


1 


An 


6 


q 


m 


+ —-1+4^-lW --mg 5-^ u“V«V^4> 


m^ 


9" 


1 


m^ 


+ r~ dr'G^^^ + f" dr'm'u^' + [" dr'G^^^ 


m 


- -mg ( - I u^u^w^'' + 




VA«h / ” dr’Gf^ys'^^'u^' 


+ r dT%'A>G^%,s'U^'u^' + r dr' - 2u^VaGf^^u'^) 


' —OO 

rr- 


+ I dr' Gf^y - qVxA> / dr' iw>^''w^^Gfl + w^^Gffu^u^ 


VI- 


Ad- 


(64) 


describing the motion of a point particle under the inhuence of the hrst-order coupled self¬ 
force in a scalarvac spacetime. Note that we have adopted the “dot” notation to distinguish 
the oh-diagonal components of Gab from the diagonal ones. Despite its complicated appear¬ 
ance, several terms carry simple interpretations. The leading order local term is simply the 
acceleration in the background spacetime. The second local term coming from 

the order-reduced acceleration and the Ricci tensor, may be interpreted as a background 
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force on a modified mass grfi. The third term is a combination of the order-reduced self¬ 
acceleration’s time derivative and the the mixing tensor The expression above differs 

slightly from the one found in [1] where a small error was made in the Fermi coordinate 
components of the gradient of the metric perturbation. Specihcally, the component Vahu 
in [1] is missing a term given by —2g<haa. We also find that the tail has a much more 
complicated structure than in the vacuum situation. The rich structure of the tail can be 
seen in the results of [1], which use the regular solutions 

0[tail] := 4 f {m — q^)G', {x, z)u^u’^ dr + q f G',{x,z)dT, (65) 


and 


h„/ 3 [tail] := 4 / (m - g<h)A z)uf"u'" dr + q A^^^G^s\-{x, z) dr, (66) 


which both contain off-diagonal contributions from the Green function Gab- 

VII. ELECTROVAC SELF-FORCE 

A. Perturbed electrovac action 

For our second example we consider a massive, charged particle moving in an electrovac 
spacetime. Here we consider a background spacetime with metric gafs and a background 
vector held with corresponding bulk background action 


Shalt = \ f 


(67) 


where Fafs is the held strength tensor of the electromagnetic held Fa/s = 2V[oA_a]. The 
stress-energy tensor of the bulk matter is given by 


= 2 - -g^gF.^F'^'' ) , 

which implies F = 0. The background held equations read 


R-aji TajS 

'^^Fa/S = 0, V[Q,F/3q = 0. 


( 68 ) 


As before the background spacetime is perturbed by the presence of the point particle 
carrying an electric charge e and mass m. The particle follows a world line in the background 
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spacetime 7 parameterized by proper time r. We denote the particle’s four velocity in the 
background spacetime by and normalize it in the background metric such that Qapu^u^ = 
— 1 . The minimally coupled world line action of the point particle is given by 

S'pp = —m j dr + e j A^u^dr. (69) 

The perturbation introduced by the particle sources a perturbation of the fields dehned as 
the difference between the full non-linear helds, Oa and ^a/ 3 , and their background values 

bfj, flp Af^ 
hajS Scip dap- 

The two perturbed fields are packaged using our multi-index notation into the doublet 
t/’A := {ba-,hap}- The variations ddr = —^hfj_^u^u''dT and 6A^ = couple to a source 
doublet ju^ = — drS^i^x, z)s^{z) given by 

:= {eti", l/2mu^u''}. 

We shall also require the hrst two derivatives of the source 

= {ea", 

and 

= {ed°‘, m } 

in the expansion of the world line. 

In order to apply the multiheld formalism of Sec. Ill, we must expand the action to 
second-order in the held perturbations and introduce a gauge. This allows us to compute 
the propagators associated with the helds 'ipA- We hnd that the second-order variation of 
the action reads 

= j (2V[„6^] Vpb^ - ^F^pF^^ 

_ - -^PF^'^^V^hVph.s + \v,Kp dV. 

(70) 

Unlike the scalarvac scenario, we are unable to hnd a gauge that eliminates the derivative 
coupling while maintaining a smooth relation to the Lorenz gauge. We therefore choose the 
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gauge to be the standard Lorenz gauge for both helds, which translates to the gauge doublet 
Ga = {1/2Lq,, Vo-a"} in the multiheld language. After simplihcation we hnd that the bulk 
action reads 

^'^buik = -^J (&« b^+\K^ sh^^yv, 

where 




^a/3jS 




yo/ 37(5 ._ g 5 ]g _ gd^p-yS 


from which we build the gauge-hxed held operator in matrix form 

-g^^n + 971“^ 0 

y{o/3)7(5y^ 1 (^J^(ag){y&) _|_ _/y^(75)M)^ _ l^a/37<5Q 


6*‘> = 


(71) 


B. Local electrovac self-force I: expansion for vanishing (§ab 


As in the scalarvac scenario we hrst establish the local self-force result in hat spacetime 
where Sab = 0. The absence of any coupling between the metric perturbation and the 
background electromagnetic held on the world line greatly simplihes the local analysis in 
comparison to the scalarvac scenario. Using Eq. (55) we hnd that the CAP variation of the 
ehective action yields 


I dr'W[aD^'^pv?' + / dr' ( - w, 

J —oo ^ J —oo 

from which we read oh the quantities 


ia-- i«Vad D’p,,,u‘‘‘u‘‘‘, (72) 


-4% = 0 . 

V"" = . 

5 = mui 

Substituting these expressions and the expressions for the time derivatives of the source into 
Eq. (43) gives the result 

47rf« = wj^dp (73) 
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At zeroth-order in S’ab, the local self-force for the electrovac spacetime is simply the sum 
of the well-known gravitational anti-damping force and the Abraham-Lorentz-Dirac force. 
The additional term that we saw in the scalarvac case is absent here because the mass of the 
world line doesn’t depend on the background electromagnetic held. Using the background 
equations of motion maa = eF^isu^ and maa = eu^u^'V + ^F^pF^'^u^, we obtain the 
order-reduced form of the equation of motion 

47rf„ = -\me (ll - 2^) u^u^V.F^p - (4-2^- 5^) F^^F^^u^ (74) 

O \ TTl j O \ 7Tl 6 / 

which is free of runaway solutions. 


C. Local electrovac self-force II: contributions from S’ab 


We now compute the local contribution to the equation of motion resulting from bulk 
curvature corrections to the Green functions. The starting point is the bulk action in 
coordinate form 


5'^buik = Sh + Sf, + S„ 


(75) 


where 

Sb = -\J dV + 2g^^ + 0(r2)), 

= ~j + 0(r2)). (76) 

Inspecting Eqs. (76) and dropping terms higher than second-order in background derivatives 
we hnd the curvature potentials to be 





gllUgal 3 




and 


^ab,i, ^ ^ 


( 4 pM]M Q 

y«/37/i _gtJ-i'^ 


\ 


1 






0 

0 
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Using the RNC expansions in Eqs. (60) and (61), the second-order constituents of S’ab take 
on the values 

(77a) 


(B» <>« + - (B“'’,1^+ByV;"'’)). 

(77b) 

O 

^oiy,y _ _ 

3 

(77c) 

^al3'y,y _ _^_pQ[7^/7]/3 _ ^a/3 p'yy 

(77d) 

o 

^al3'YS,y _ ° yuj^a^ri-y jA) 

3 ' {vT)y' 

(77e) 


(77f) 

— J^{al3){y5) 

(77g) 


where again we adopt the shorthand notation Ay := Aay°‘ for contractions with Riemann 
normal coordinates. Again we proceed by replacing powers of the RNCs with derivatives 
with respect to momentum to write ^ in the form of (45). From there, we directly proceed 
to the expression for^ in terms of the integrals in Eq. (46). 

We begin the integrations by considering the terms. We hnd that the kinetic tensor 
contracts with the integral in such a way that it vanishes geometrically due to the 

symmetries of the Riemann tensor. The result of integrating the single derivative angu¬ 
lar momentum tensor is ^Rayu^, which combines with the algebraic mass term given by 
^Rayu^, to give ^e^Rapu^ = u^, showing a deviation from the test 

field electromagnetic local self-force ^e^Ray,- The non-zero terms are given by a com¬ 
bination of the RNC expansion of the metric determinant in and the terms in 

yielding ^rn?Rayu^ and 2Faj3F^^u^, respectively. The local contribution from the 
mixing is simply AemF^/sa^. Adding the three parts gives the local electrovac self-force from 
bulk curvature insertions 

47r/„ = (^RyyF' + 2Fy„F’'pU^^ + ^e^Ryuu" + AmeFy^a’"^ . (78) 
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D. Electrovac equation of motion 


Combining the local results with the tail terms coming from the variation in Eq. (55), 
and writing in terms of the held strength tensor, we arrive at the covariant equation 


1 


1 


dvr 




1 


ma^ =eF^^''uu + - -me < 11 - 2 —^ u — -e ' 4 - 2— 5 - - 6 - 5 - 


.m? 


9 9 


+ — ( w^^^^Vx + -— FP^ 


1 e 


2 m 




dT'Gl%^,u^'u^' 


+ 


2e^ux f + emw^^uxu^ f 

J—00 J—i 


ylcr^Ajret c 


+ emwP^ux [ — em < ) [ drGyft^^iU^' (79) 


V5|/3' 


describing the motion of a charged particle in an electrovac environment. 


VIII. CONCLUSION 

In this paper we have derived the hrst-order equations of motion for a small compact 
body moving in various non-vacuum spacetimes using techniques from effective held theory. 
We began with the simple example of a scalar particle with N charges coupled to an N- 
component scalar held, which allowed us to introduce the ehective action procedure. With 
the toy model we introduced several technical methods such as the local expansions and 
dimensional regularization that were used throughout the paper to obtain the force at the 
location of the particle, while keeping within the simplest possible context. We also intro¬ 
duced our procedure for building the curved spacetime propagator as a series constructed 
from the propagator in hat spacetime and vertex corrections. Using analogies with the toy- 
scalar model, we then derived a general expression for the self-force for general non-vacuum 
spacetimes using multiheld notation. With this notation, we gave an explicit expression for 
the local self-force in terms of quantities that can be read oh directly from the action along 
with a set of integrals falling under a single master integral expression. 

We then applied the formalism to derive equations of motion for point particles in scalar- 
vac and electrovac spacetimes as the main results. In the scalarvac spacetime, the mass 
of the world line evolves with the background scalar held, which led to new local world 
line couplings in the equation of motion relative to the decoupled scenario. In contrast, 
the mass of the world line in the electrovac case is conserved and we found that the local 
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self-force at zeroth order in the bulk vertex S’ab was simply the sum of the vacuum and 
test electrovac results. The coupling between held perturbations in the bulk does create 
additional world line terms, however, and the local force in the presence of curvature cannot 
be written as a simple sum of the decoupled local force terms. In both the scalarvac and 
electrovac cases the non-local tail terms were written as integrals involving the gravitational 
and scalar/electromagnetic propagators and additional non-diagonal Green functions. These 
non-diagonal Green functions contribute both local terms and tail terms. We obtained their 
local contributions by building an approximation using the hat-spacetime diagonal Green 
function with a series of mixing insertions and found that a single mixing insertion was 
satisfactory to obtain the local force. 

In this paper we limited our focus to the ehective action for the world line and left the 
computation of the regular held for future investigations. This has the beneht of ohering 
a direct route to the world line without the need of solving the held equations. Instead, 
we use easy-to-work-with Green functions dehned in hat spacetime and treat curvature 
corrections perturbatively as interaction vertices. The drawback is that one must do a 
separate computation of the ehective action with an external held line to obtain the retarded 
held. The beneht of using the held equations is that the solutions for the helds are a 
byproduct of the self-force calculation instead of a separate endeavor. With the solutions 
themselves, at least locally, one can then use the Detweiler-Whiting routine to dehne a 
singular held from which one can derive regularization parameters for a practical calculation. 
The Detweiler-Whiting routine has been applied to an EFT calculation of third-order scalar 
perturbations, but results for regular metric perturbations are still lacking. 

The equation of motion results presented here and in [1], in combination with the analysis 
in [2], promise to open new avenues of research in the self-force program. One future direction 
is exploring the self-force equations of motion in alternative theories of gravity. Fundamental 
scalar helds are predicted to exist in many low energy limits of quantum gravity theories. 
Such scalar-tensor theories can be written in the form of (47) with an appropriate choice 
of conformal frame. We are currently formulating the self-force problem for general small 
bodies in these theories using methods present here as well as those in [1] with eventual 
plans to explore the weak-held limit. Additionally, we wish to revisit Hubeny’s overcharging 
scenario. Using the coupled electrovac equation of motion, one could expand the local 
analysis conducted by Hubeny, which ignored the tail piece, to include the new local terms 
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generated by the coupling as a preliminary condition for potential overcharging candidates. 
The candidates could then be followed up with a full numerical self-force calculation via the 
regular field using the results of [1] and [2], 
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Appendix A: Fixing the gange 


Under inhnitesimal gauge transformations the metric perturbation undergoes a change 

Let us now check whether the modihed Lorenz gauge condition for the scalarvac problem 
Ga = La + a<h;o/ for constant a alters the path integral in any non-trivial way. We consider 
the gauge transformation 

where 

= L>^{h) + Dx^ + 

The functional determinant of G°‘{h^)/Sx^ = is independent of the dynamical 

helds hap and / and therefore contributes nothing to correlation functions. Now, as the 
determinant is independent of hap and /, we can remove it from the path integral and 
relabel the integration variable such that 

j VhVfe^^ = det + j '^X j DhDfG^^^L]s^G{h,f)). 

Generalizing the gauge condition to read 

G^(h,/)^G^(h,/)-A^(a;) 


we hud 


VhVfR^ = det (r^n + R“^) Vx + A^). 
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Since the above equation holds for arbitrary functions X^{x), it must also hold for any 
normalized linear combination of such functions. We choose the normalization N{^) with 
respect to a Gaussian of width ^ 


VhVfe^^ = det iV(0 J VXe~^^ J Vx j + a/$’^ - A^) 

= ( det + i?“^) iV(0 [ 


choosing ^ = k and we see that we have recovered the gauge fixing term in (27) 


[1] P. Zimmerman and E. Poisson, Gravitational self-force in nonvacuum spacetimes, Phys. Rev. 
D 90, no. 8, 084030 (2014), arXiv;gr-qc/1406.5111. 

[2] T. M. Linz, J. L. Friedman and A. G. Wiseman, Combined gravitational and electromagnetic 
self-force on charged particles in electrovac spacetimes, Phys. Rev. D 90, no. 8, 084031 (2014) 
arXiv;gr-qc/1406.5112. 

[3] Y. Mino, M. Sasaki, and T. Tanaka, Gravitational radiation reaction to a particle motion, 
Phys. Rev. D 55, 3457 (1997), arXiv;gr-qc/9606018. 

[4] T. C. Quinn and R. M. Wald, Axiomatic approach to electromagnetic and gravitational ra¬ 
diation reaction of particles in curved spacetime, Phys. Rev. D 56, 3381 (1997), arXiv:gr- 
qc/9610053. 

[5] S. E. Gralla and R. M. Wald, A rigorous derivation of gravitational self-force. Glass. Quantum 
Grav. 25, 205009 (2008), arXiv:0806.3293. 

[6] A. Pound, Self-consistent gravitational self-force, Phys. Rev. D 81, 024023 (2010), 
arXiv;0907.5197. 

[7] L. Barack and N. Sago, Gravitational self-force on a particle in eccentric orbit around a 
Schwarzschild black hole, Phys. Rev. D 81, 084021 (2010), arXiv; 1002.2386. 

[8] L. Barack and N. Sago, Beyond the geodesic approximation: conservative effects of the grav¬ 
itational self-force in eccentric orbits around a Schwarzschild black hole, Phys. Rev. D 83, 
084023 (2011), arXiv:1101.3331. 

[9] N. Warburton, S. Akcay, L. Barack, J. R. Gair, and N. Sago, Evolution of inspiral orbits 
around a Schwarzschild black hole, Phys. Rev. D 85, 061501 (2012), arXiv:1111.6908. 


38 



[10] N. Warburton and L. Barack, Self force on a scalar charge in Kerr spacetime: circular equa¬ 
torial orbits,, Phys. Rev. D 81, 084039 (2010) 

[11] N. Warburton and L. Barack, em Self force on a scalar charge in Kerr spacetime: eccentric 
equatorial orbits, Phys. Rev. D 83, 124038 (2011) [arXiv;1103.0287 [gr-qc]]. 

[12] S. Detweiler, Gravitational radiation reaction and second-order perturbation theory, Phys. Rev. 
D 85, 044048 (2012), arXiv:1107.2098. 

[13] S. E. Gralla, Second-order gravitational self-force, Phys. Rev. D 85, 124011 (2012), 
arXiv; 1203.3189. 

[14] A. Pound, Second-order gravitational self-force, Phys. Rev. Lett. 109, 051101 (2012), 
arXiv;1201.5089. 

[15] A. Pound, Nonlinear gravitational self-force: Field outside a small body, Phys. Rev. D 86, 
084019 (2012), arXiv:1206.6538. 

[16] A. Pound and J. Miller, A practical, covariant puncture for second-order self-force calculations 
(2014), arXiv;1403.1843. 

[17] P. Diener, I. Vega, B. Warded and S. Detweiler, Self-consistent orbital evolution of a particle 
around a Schwarzschild black hole Phys. Rev. Lett. 108, 191102 (2012) arXiv;1112.4821. 

[18] B. Warded, C. R. Galley, A. Zengino?lu, M. Casals, S. R. Dolan and A. C. Ottewid, Phys. 
Rev. D 89, no. 8, 084021 (2014) [arXiv; 1401.1506 [gr-qc]]. 

[19] A. Le Tiec, A. Buonanno, A. H. Mroue, H. P. Pfeiffer, D. A. Hemberger, G. Lovelace, L. E. 
Kidder, M. A. Scheel, B. Szilagyi, N. W. Taylor, and S. A. Teukolsky, Periastron advance in 
spinning black hole binaries: Gravitational self-force from numerical relativity, Phys. Rev. D 
88, 124027 (2013), arXiv:1309.0541. 

[20] A. G. Shah, J. L. Friedman, and B. F. Whiting, Finding high-order analytic post-Newtonian 
parameters from a high-precision numerical self-force calculation, Phys. Rev. D 89, 064042 
(2014), arXiv;1312.1952. 

[21] E. Poisson, A. Pound, and 1. Vega, The motion of point particles in curved spacetime. Living 
Rev. Rel. 14, 7 (2011), arXiv:1102.0529. 

[22] V. E. Hubeny, Overcharging a black hole and cosmic censorship, Phys. Rev. D 59, 064013 
(1999), arXiv;gr-qc/9808043. 

[23] T. Jacobson and T. P. Sotiriou, Overspinning a black hole with a test body, Phys. Rev. Lett. 
103, 141101 (2009), arXiv;0907.4146. 


39 



[24] E. Barausse, V. Cardoso, and G. Khanna, Test bodies and naked singularities: Is the self-force 
the cosmic censor?, Phys. Rev. Lett. 105, 261102 (2010), arXiv:1008.5159. 

[25] M. Colleoni and L. Barack, Overspinning a Kerr black hole: the effect of self-force, 
arXiv;1501.07330 [gr-qc]. 

[26] P. Zimmerman, I. Vega, E. Poisson, and R. Haas, Self-force as a cosmic censor, Phys. Rev. 
D 87, 041501 (2013), arXiv:1211.3889. 

[27] S. E. Gralla, Motion of small bodies in classical field theory, Phys. Rev. D 81, 084060 (2010), 
arXiv; 1002.5045. 

[28] S. E. Gralla, Mass, charge, and motion in covariant gravity theories, Phys. Rev. D 87, 104020 
(2013), arXiv; 1303.0269. 

[29] T. G. Quinn, Axiomatic approach to radiation reaction of scalar point particles in curved 
space-time, Phys. Rev. D 62, 064029 (2000) gr-qc/0005030. 

[30] S. Detweiler and B. E. Whiting, Self-force via a Green’s function decomposition, Phys. Rev. 
D 67, 024025 (2003), arXiv:gr-qc/0202086. 

[31] W. D. Goldberger and I. Z. Rothstein, An Effective field theory of gravity for extended objects, 
Phys. Rev. D 73, 104029 (2006) arXiv;hep-th/0409156. 

[32] G. R. Galley, Classical Mechanics of Nonconservative Systems, Phys. Rev. Lett. 110, no. 17, 
174301 (2013) arXiv; 1210.2745. 

[33] G. R. Galley, D. Tsang and L. G. Stein, The principle of stationary nonconservative action 
for classical mechanics and field theories, arXiv; 1412.3082 [math-ph]. 

[34] G. R. Galley, B. L. Hu, and S.-Y. Lin, Electromagnetic and gravitational self-force on a rela¬ 
tivistic particle from quantum fields in curved space, Phys. Rev. D 74, 024017 (2006), arXiv;gr- 
qc/0603099. 

[35] G. R. Galley and B. L. Hu, Self-force on extreme mass ratio inspirals via curved spacetime 
effective field theory, Phys. Rev. D 79, 064002 (2009), arXiv;0801.0900. 

[36] G. R. Galley, A nonlinear scalar model of extreme mass ratio inspirals in effective field theory 
I. Self force through third order. Glass. Quant. Grav. 29, 015010 (2012) arXiv; 1012.4488. 

[37] G. R. Galley, A Nonlinear scalar model of extreme mass ratio inspirals in effective field the¬ 
ory II. Scalar perturbations and a master source. Glass. Quant. Grav. 29, 015011 (2012) 
arXiv; 1107.0766. 


40 



[38] C. R. Galley and M. Tiglio, Radiation reaction and gravitational waves in the effective field 
theory approach, Phys. Rev. D 79, 124027 (2009) arXiv:0903.1122. 

[39] J. Schwinger, Brownian Motion of a Quantum Oscillator, J. Math. Phys. 2, 407 (1961). 


41 



